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Abstract. The identification of dark matter in our particle physics model is still a 
very open question. Here we will argue that axinos can be successful dark matter 
candidates in models with supersymmetry and the axion solution of the strong CP 
problem. Axinos can be the lightest supersymmetric particle (LSP), or can be heavier 
than the LSP. Axinos can be produced in the right abundance by thermal scatterings 
and if they are the LSP also by out of equilibrium decays of the lightest superpartner of 
SM fields (LSPSMs). On the other hand heavier (not LSP) axinos can generate a part 
of the neutralino LSP dark matter. Depending on the nature of the supersymmetric 
spectrum, and if R-parity is strictly conserved or slightly broken, very different signals 
of the LSP axino scenario can arise at colliders and in astrophysics. 
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1. Introduction 

The nature of dark matter (DM) as a particle is still unknown today, since its main 
evidence relies only on the gravitational interaction and is universal. On the other 
hand, from some of the present data and numerical simulations of structure formation, 
we do know that it must be a neutral, cold, very probably collisionless (i.e. quite 
weakly interacting) and very long lived particle [1] . Unfortunately a particle with these 
characteristics is not contained in the standard model (SM) of particle physics: the only 
neutral stable and massive candidates, the electroweakly active neutrinos, are so light 
that they are at most hot DM and therefore only a subdominant component. DM has 
therefore to be part of a larger picture and of any physics beyond the SM. 

The probably best motivated models of this kind rely on supersymmetry (SUSY), 
which is the unique extension of the Poincare' algebra and calls for a doubling of all 
degrees of freedom with spin difference As = ±| [2]. In this context it is then clear 
that more particles can be suitable DM candidates, if they are the lightest one and 
sufficiently long-lived, in particular the very well studied cases of the neutralino or the 
gravitino. 

But if we invoke the Peccei-Quinn (PQ) solution to the strong CP problem in 
super symmetric models, a new multiplet has to be introduced, the axion multiplet [3lll]. 
Such a multiplet must by its nature interacts with the SM particles, but the scale 
of its interaction is suppressed by the scale at which the PQ symmetry is broken. 
Fa- Therefore, the fermionic component of the multiplet is naturally a very weakly 
interacting particle and can easily be the lightest state of the spectrum, but it can also 
be heavy. We will present in this paper a summary of the axino cold DM (CDM) 
scenario [5], El [TJ El [HI [ini [H] and explore the implications of axino DM for the 
phenomenology of supersymmetric models and for future indirect detection of DM. 

2. Axino models and axino mass 

In discussing axino models, one should refer to the corresponding axion models. So, let 
us start with the axion shift symmetry and the reparametrization invariance as discussed 
in [12]. 

The PQ solution of the strong CP problem requires the introduction of the axion a, 
which renders the 6 parameter dynamical and allows it to relax to zero after the QCD 
phase transition. An axion a is a pseudoscalar boson coupling to the gluon anomaly as 

^0 = ^ G,.G''^ , (1) 

where the dual field strength is G'^'^ = ^e'^'^^^Gpa^ without any other interaction term in 
the potential V. Below the QCD chiral symmetry breaking scale an axion potential is 
developed, which arises purely from integrating out the strongly interacting fields with 
that anomaly term. If the original potential contains other axion dependent terms, they 

I Below, 9 denotes the conventional 9 = 9o + Arg.Det.mq. 
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should be extremely small and not affect the position of the minimum, such that the 
axion v.e.v. \6\ ~ (a) /Fa should be extremely small, < 10"^^. 

Let us now focus on the QCD interactions containing the axionic degree of freedom. 
The axion effective Lagrangian below the chiral symmetry breaking scale reads 

- (qL m qne'^^' + h.c.) + Cs^G^l^G"^'^ (or ^dct) (2) 
9 

where 9 = a/fs with the the axion decay constant fs defined up to the domain wall 
number {fs = NdwFo) and q are the SU(3)c charge carrying quark fields. The ci term 
is the derivative coupling with quarks respecting the PQ shift symmetry, the C2 term 
is related to the phase in the quark mass matrix, and the C3 term is the anomalous 
coupling or the determinental interaction £dct- -Cicptons.e is the axion interaction with 
leptons, which in principle can contain other constants c\,c\. The coupling constants 
ci,C2, and C3 are obtained below the axion scale fs after integrating out the heavy 
degrees of freedom responsible of the PQ symmetry breaking. The mass parameter m 
is defined to be real and positive below the electroweak scale. 

The Lagrangian ([2]) has a shift symmetry a a+ (constant), which reparametrizes 
the couplings Ci,C2, and C3. Explicitly, changing the phases of the quark fields 
qi — > e'"°^^''gL and q^ e~*""^^''g/?, we obtain the following reparametrization, where 
the effective one point irreducible action Vipj[a{x), A^(x); Ci, C2, C3, m, Aqcd] changes to 

Tipi[a{x),A"^{x);ci - a,C2 - 2a, C3 + 2a, m, Aqcd]- (3) 

So we see immediately from this transformation for a single quark, that if it is massless, 
the corresponding Ci, C2 parameters disappear and we can shift away the anomaly term 
completely with no physical effect. This is in fact one alternative solution to the strong 
CP problem, and see |I2] for a detailed discussion. 

For determining the axion mass, all ci, C2 and C3 terms may be relevant, but only the 
combination C2 + C3 actually appears [12]. Usually, in the field theoretic axion models, 
we start with ci = 0. In any case, note that the ci term can be reabsorbed in the C2 
term using integration by part and the quarks equations of motion. So, in the next 
sections we just start with the couplings C2 and C3. 

Usually, Fa is defined by transferring all couplings of the axion to the coefficient 
of GG and rescaling C3 to one. On the other hand, fs is defined to be the VEV of the 
singlet field a breaking the PQ symmetry. It turns out that C2 + C3 is an integer, not 
necessarily one in the pseudoscalar field space and it determines the number of minima 
in the axion periodic potential. Thus, this integer is called the domain wall number 

ndw m 

NdW = |c2 + C3I = Tr QpQ (^/'colored)^(V'colored) , (4) 

where the trace is taken over all colored fermions 'i/'coiored, ^ is the index of their SU(3)c 
representation and the PQ charge Qpq is given for the left-handed chiral representations. 
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The scale Fa is constrained by astrophysical and cosmological bounds to lie in the 
narrow axion window 10^°GeV < < lO^^GeV [T2] . 

Note that above the electroweak (EW) symmetry breaking scale in principle also 
couplings to the EW gauge bosons and the Higgs fields may arise and we have to write 
then the effective Lagrangian as 



Then we can define an extended shift symmetry, including also transformations of 
the EW charged fields; these can be changed independently to the colored degrees of 
freedom, so to have c^^ew = and leave only the anomalous coupling to the hypercharge 
gauge bosons. As in the case of QCD, such a coupling could be shifted away completely 
if one of the leptons were massless; since the electron mass is quite small, the residual 
effects, contained in the c{ 2 terms, is negligible for many practical purposes. 

2.1. Axion Models 

There are several types of C2 and C3 couplings which define different axion models. If 
C2 = and C3 7^ due to the existence of PQ charge carrying heavy quarks, the model 
is called the Kim-Shifman-Vainshtein-Zakharov (KSVZ) model. If the coupling C2 is 
provided by the electroweak scale Higgs doublets, while C3 = 0, it is the Peccei-Quinn- 
Weinberg-Wilczek (PQWW) model. If the phase C2 is provided by an electroweak singlet 
with C3 = 0, it is the Dine-Fischler-Srednicki-Zhitnitski (DFSZ) model. The model- 
independent axion in superstring models give instead C2 = and C3 = 1. These values 
enable us to write down the axion-nucleon-nucleon couplings unambiguously for each 
model [l2]. But in general axion models may contain both C2 and C3 with C2 + C3 7^ [I2] 
and may have the family dependencies of the variant axion [H] or invisible axion |15j . 

2.2. Axino, SUSY breaking and axino mass 

In the case of a supersymmetric model, the axion field is the pseudoscalar part of a 
whole chiral multiplet $. Note, however, that the reparametrization invariance Eq. ([2]) 
still holds and represents a freedom in choosing the Ci, C2, and C3 terms. We choose here 
the basis where the C2 term is transferred to the C3 term, and hence $ interaction is 



where now $a = (s + ia) / \/2 + -i^a + (F term) is the chiral multiplet containing 
the saxion s and axion a and their fermionic partner the axino a, while Wq, is the 



i^e,eff>EW 




(5) 




(6) 
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vector multiplet containing the gluino and the gluon field strength, and W"Wq, l^^i = 
-2zAV"'a;„A"-iG;i^G'"'^'^ + iG';^^G'"^^ + D2. Here, a, is the QCD couphng constant. An 
analogous interaction is present for the hypercharge gauge multiplet with the additional 
coupling C3,y as discussed previously. 

As long as SUSY is unbroken, the axion multiplet remains light, since it is protected 
by the U(1)pq symmetry fi6[ [T71 [T8] . This symmetry implies that no supersymmetric 
mass parameter is allowed for the axion multiplet since, as discussed above, the axion 
does not have a potential V (i.e. terms in the superpotential W with SUSY). 

Both saxion and axino masses are split from the almost vanishing axion mass if 
SUSY is broken, either at tree level via the v.e.v. of some scalar field in the model 
and mixing with the other neutralinos or via loop diagrams involving multiplets with 
split masses. The precise value of the axino mass depends on the model, specified by 
the SUSY breaking sector and the mediation sector to the axion supermultiplet. Most 
probably, the saxion mass is around the soft mass scale MgusY- The axino mass should 
also be near this scale as well. But the axino mass can also be much smaller than 
^susY [H [13 [H] or much larger than Msusy [IE]- Therefore, we take the axino mass 
as a free parameter. 

If R-parity is not conserved, the lightest supersymmetric partner of the SM particles 
can decay to ordinary particles. If R-parity is conserved, it cannot decay to ordinary 
SM particles, but it can decay to axino or/and gravitino if they are lighter. Thus, the 
axino cosmology depends crucially on the R-parity realization. Here, we consider first 
models with R-parity conservation and the thermal history of the universe can be very 
different depending on the hierarchy between the axino mass and the mass of the LSP 
(of SM multiplets) M^. Firstly, we consider the case < and next > M^. 

The cosmology of a weakly-interacting massive particle (WIMP) and an extra- 
WIMP depends on several temperatures. For example, the neutralino cosmology 
depends on the neutralino freeze-out temperature fi9[ [20] and the gravitino/axino 
cosmology on the reheating temperature after infiation [211 [22] . We therefore define 
the following temperatures relevant for the axino cosmology: 

Ta-dcp = axino decoupling temperature 
Tr = reheating temperature after infiation 

Tfr = neutralino freeze — out temperature (7) 
Ta=rad = axiuo — radiatiou equality temperature 
Td = radiation temperature at x or a decay 

where note that Td corresponds to a different temperature for rria < > M^. 

3. Axino cosmology with < 

Let us consider the axion supermultiplet together with the Minimal Supersymmetric 
Standard Model (MSSM) fields. Then the lightest supersymmetric particle in the 
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MSSM (LSPSM) X plays an important role. The decoupling temperature of the axino 
supermultiplet is of the order [23] , 

Cosmology with the saxion s is a simple extension of the standard cosmology if the 
saxion mass is around the SUSY breaking scale [2lj or larger [25], but its effect is not 
so dramatic as the effect of the axino. If the axion interaction ever was in thermal 
equilibrium, e.g. Tr > Ta-dcp, a substantial axino number density survives to the 
present day and axinos have to be very light. Axinos with mass in the eV range from 
this epoch have been considered as hot DM [16] or warm DM for masses in the keV 
range [23]. In the gauge mediated SUSY breaking (GMSB) scenario, the gravitino is 
probably the LSP and the possibility of primordial axinos decaying to gravitinos has 
been considered as well |26]. Here let us focus on the CDM axino LSP scenario. 



3.1. Producing axinos in the Early Universe 

We briefly review here the two main mechanisms that produce axinos in the early 
universe. In principle also other sources could be present like Q-balls decay |27j . 
We concentrate here on the hadronic type of axion models and define the axion 
supermultiplet in the basis where the C2 term is zero. Then the axino does not 
interact directly with the MSSM multiplets apart from the gluon and hypercharge 
vector multiplets and does not mix substantially with the standard neutralinos or other 
fermions. So we neglect any interaction with the leptons or the EW gauge bosons, that 
may appear in the DFSZ type of models, and that can only increase the production 
cross-section. It is worthwhile to recall here again that the shift between C2, C2 and 
C3, C3,EW couplings is simply a matter of definition of the axion interactions as long as 
it is nearly a mass eigenstate and if we choose C2,C3^ew = 0, the other coupling is 
suppressed by the leptons Yukawa and therefore negligible for axino production. 



3.1.1. Thermal scatterings Any particle, even very weakly coupled, is produced in the 
thermal plasma by scatterings of the particles that are in thermal equilibrium. As we 
have seen axinos couple directly to the gluons and gluinos via the "anomaly" coupling 
in Eq. ([6|), i.e. in components 

i^-a,, = ^^^s^'^'-A^G;:, , (9) 

where A'' is the gluino field. So many scatterings in the primordial plasma involving 
colored particles produce axinos The axino number density is given by solving a 
Boltzmann equation of the type 
dn~ 

+ ?,Hna = Y^{a{i + j ^ a + . . .)vrei)ninj + ^(r(i ^ a + . . .))ni , (10) 



§ The same happens also in the case of the gravitino, but with different vertex structure and scale |28j . 
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Figure 1. Maximal reheating temperature as a function of the axino mass obtained 
by requiring that the axino energy density is below the present DM density [5]. The 
difference between solid and dashed lines is due to the inclusion of the decay term in 
the Boltzmann equation (jlOp . In the yellow area we expect the non-thermal production 
via out of equilibrium decays to be also substantial. 



where we are neglecting back-reactions, that are suppressed by ria -C Ui. At high 
temperature the 2-body scatterings dominate the r.h.s., since they contain a vertex 
given by the dimension 5 operator in Eq. ([6]) and show a characteristic hnear dependence 
on T. So most of the axinos are produced at the highest temperature, and the axino 
number density is proportional to that temperature, which we take to be Tr. Some of 
the two body scatterings are IR divergent due to the massless gluon propagator; in the 
thermal bath such a divergence is screened by the presence of a thermal gluon mass 
~ gT. In our computation we introduced such IR cut-off by hand [6] . A self-consistent 
procedure is instead to perform a full resummation of the hard thermal loops as done 
in [29]. In general we expect 0{1) corrections from higher orders terms in as, especially 
at low temperature [29]. There as well the decay terms start dominating and the number 
density is no more proportional to the reheating temperature, but depends instead on 
the supersymmetric spectrum, in particular the gluino and squark masses [7] . Using the 
expression for the present axino energy density as 

ma. — 77:;;^ = 0.72eV , 11 

s(T) V 0-2 / 

where s{T) = 2.89 x 10^ ( „ ^Li^- ) cm~'^ is the present entropy density, we can then obtain 



,2.726X 

a bound on the reheating temperature in Fig. [H 
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3.1.2. Out of equilibrium decays An axino population is also generated by the LSPSM[iJ] 
decay after it freezes out from the thermal bath. The heavier superpartners cascade- 
decay quickly into the LSPSM (or very rarely to the axino LSP itself as we discussed 
above) while still in equilibrium, but the LSPSM has a lifetime much longer than its 
freeze-out time: since the axino couplings are suppressed by the PQ scale Fa — 10^^ 
GeV, the LSPSM lifetime for 100 GeV mass is of the order of seconds. Then the freeze- 
out process is unaffected since the decay takes place only much later. 

In this case, thanks to R-parity conservation, the axino energy density can be 
directly computed from the LSPSM would-be-relic density as 

^5.'^ = ^LSPSM- (12) 

^LSPSM 

If the mass ratio is not too small, we still have a connection with the classical WIMP 
mechanism in case the LSPSM is a neutralino. On the other hand in this scenario the 
LSPSM can be more generally any superpartner which may freeze-out with a sufficiently 
large number density. 

A couple of problems can arise if the LSPSM decay happens too late: 

• Big Bang Nucleosynthesis (BBN) can be spoiled by the energetic "active" particles 
produced in the decay along with the axino: the strong limits on the injection of 
energetic particles depend on the electromagnetic/hadronic nature of the produced 
showers, the LSPSM number density and its decay time [22]. In general such limits 
are weak for the axino case since the LSPSM lifetime (excluding a strong mass 
degeneracy) is below lO^s, but they can affect the region of small mass for both 
the neutralino and stau LSPSM [HllTj. This constraints disappear easily for colored 
LSPSM like the stop since the number density is reduced also by the Sommerfeld 
enhancement and the lifetime is very short as long as the decay to top is allowed |llj . 
Therefore a stop LSPSM is perfectly viable, but a very large stop mass and axino 
mass of a few TeV is needed to produce the whole DM density. 

• Are axino from the decay cold enough to be CDM? They may be relativistic at 
production even if the LSPSM is not and they in general have a non-thermal 
spectrum. Their velocity can be estimated as: 

v{T) = P^^ "^^^^ ( ^, (13) 
m-a 2m~a \g^{TB) J Td 

where Td here is the temperature of the LSPSM decay time. Axino must therefore 
have sufficient time to cool down before structure formation begins. In [31] such 
constraints have been studied and the conclusion is that an axino mass of at least 
1 GeV is probably needed if the whole DM population is produced by out of 
equilibrium decay of a LSPSM of a 100 GeV mass. 

II In passing, we point out that the gravitino problem |21| is absent if the gravitino is the NLSP, 
rua < 12 < , since a thermally produced gravitino would decay into an axino and an axion which 
do not affect the BBN produced light elements [30| . 
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Depending on the parameters and Tr, either production mechanism can give 
sufficient axinos to explain the present DM density. Once more information about 
the SUSY spectrum is available from LHC it may be possible to determine which 
contribution dominates and restrict the range of Tr [32]. Of course another possibility 
is that the axino is so light to be a subdominant (warm or hot) DM component. In the 
last case in our scenario the axion [12] could be the DM. 

4. Axino cosmology with > 

Now, let us consider the axino mass region, > M^. Here, we are interested in the 
case where the CDM density is determined by the axino and in particular the axino 
energy density is dominating the evolution history of the Universe. This is possible not 
only in the near past if the axino has not decayed yet and is DM as discussed in Sec. El 
but also if a heavy axino decayed into the DM at an earlier epoch as will be considered 
below. 

Also for heavy axino, the axino density before decay can be estimated from Ta_dcp 
or Tr as discussed above. Even in the second case, when the axion coupling never was 
in equilibrium, an early cold axino DM domination may have happened if the number 
density was sufficiently large, i.e. if Tr was larger than defined by the equality of 
axino and radiation energy density at decay: 

Im.Y.iTD = Td . (14) 

So for any Tr > T™° axinos dominate the evolution of the universe before they 
decay and produce a non-negligible amount of entropy diluting the existing number 
densities. We recall here that in SUSY theories we must always consider a relatively 
small reheating temperature 10^"^ GeV due to the gravitino problem [211 [22]. The heavy 
axino cosmology must also satisfy this upper bound on the reheating temperature. 

Heavy axinos cannot be the LSP and decay to the LSP plus light SM particles. 
This possibility was considered briefly in studying cosmological effects of the saxion in 
[25], and a more complete cosmological analysis has been presented in [10] which will 
be sketched here. Here the axino or the axino-decay produced neutralino is supposed 
to constitute the CDM fraction f2cDM — 0.23 of the universe. 

For a heavy axino decaying to a neutralino, we present a Tr vs. nia plot for 
Fa = 10^^ GeV in Fig. [2j For other parameters, we refer to [10]. The region Tr > Ta-dcp 
is above the dashed blue line. The axino lifetime is greater than 0.1 s in the red shaded 
region in the LHS and there BBN constraints may apply. The blue shaded region in 
the RHS is where the axino decays before the neutralino decouples (Tq > Tfr). The 
magenta lines (horizontal) are the contours of the entropy increase due to the axino 
decay, r = Sf/So. Above the r = 1 line axinos dominate the universe before they decay. 
The green lines (vertical) show the values of {(TannVrei), where aann is the neutralino 
annihilation cross section, in units of GeV~^ which are needed to give the right amount 
of neutralino relic density after Tfr- In Fig. [21 we use neutralino and gluino masses as 
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Figure 2. The Tr vs. ma plot for = 100 GcV and Fa = 10^^ GcV 



100 GeV and rUg = 2 TeV, respectively. For a larger Fa and a heavier neutralino 



mass, the green lines move to the right [10]. 



5. Axino DM and R-parity breaking 

The axino is a good DM candidate even if R-parity is not exactly conserved. In 
fact its lifetime can be very long, thanks to the PQ scale suppression and a small 
R-parity breaking. Axino CDM with R-parity breaking has been considered in [33] and 
subsequently in the context of the Integral anomaly in [511 ES] • Different decay channels 
and lifetimes are possible for the axino CDM depending on the R-parity breaking model 
and the axion model. If we consider bilinear R-parity breaking of the form [36] 

Wjjt^ = e.fxLiHu , (15) 

where fi is the Higgs bilinear term, and the corresponding SUSY breaking soft term 

C-softijL = BifiLiHu , (16) 

and restrict ourselves to hadronic axion models, the axino decay arises from the anomaly 
coupling with the photon vector multiplet and the neutralino-neutrino mixing generated 
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by the sneutrino v.e.v (z/j). The decay rate then reads 

2(4vr)3F2 

where = {i'i)/v with the Higgs w.e.w. given by = 174 GeV and 



■p _ '-^a77"em'"a ^2|rr |2 



f^7Z = ^^E 



7^ " m 



is the photino-Zino mixing parameter, containing the neutrahno mixing matrix S and 
the mass eigenvalues m^^ as in the case of the gravitino [37]. The hfetime is then given 
by 

For larger axino masses also the decay into Z bosons opens up, and quickly dominates 
since it is proportional to the U^^ mixing, which can be as large as one. The decay into 
W instead does not appear since we have shifted away the SU(2) anomaly vertex. Note 
that the constraints from the diffuse gamma-ray background require a very small value 
for the parameters and therefore a very small R-parity breaking, smaller than in the 
case of the gravitino DM [38] • In fact the EGRET diffuse flux already limits the lifetime 
of a DM particle decaying into a single gamma line to be larger than r > 7 x lO^^s [39], 
between 0.1 -10 GeV, with very weak dependence on the mass scale, while for lower 
masses bounds of the order of lO^^s are obtained from X-ray data [ID]. 

If the axion model is instead of the DFSZ type, direct couplings with the leptons 
and Higgses arise from the /i or ejyU terms. Then a direct mixing between axino-neutrino 
and axino- higgsino appears and other channels open up, in particular the 3-body decay 
into neutrino and via intermediate Z, which could contribute to the electron fiux|5]|. 
The mixing of the axino with leptons may arise actually also in KSVZ models at the 
one loop level [8]. Taking the mixing between the axino and Higgs multiplet typically 
as v/Fa and the Higgsino mixing to neutrino from the sneutrino v.e.v. as above, we 
obtain for this channel 



a 



(10 



26, 



if Fa \-Vlf^Hzle 



10 GeV/ VlO^i GeV/ \ 10 



-12 



where is the Fermi constant and U^j^ the mixing between Higgsino and Zino 
neutralino defined in an analogous way as U-^- The branching ratios in the different 
leptons and quarks are determined by the axino mass and the Z couplings. In general the 
same decay can also arise via the R-parity breaking leptonic trilinear coupling XLLE'^, 

% Note that also the diagram with ofF-shell photon spHtting into an electron position pair generate 
this channel, but it is then subleading in comparison to the ji' two-body channel due to aem and 
phase-space suppression. 
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and then the decay rate is given by 



24(87r)2F2mi 

-,25 



2 

R 



(21) 



^ ^ UOGeVy UOOGeVy UQiiGeV 



10 



-12 



in the hmit mg^ ^ m^, where we have taken (erug^/Fa as the effective couphng of the 
axino to the electron multiplet. Which of the two leptonic diagrams dominate depends 
on the R-parity breaking and axion model parameters. 

Recently, the PAMELA satellite experiment reported a significant positron 
excess [H], but no antiproton excess [l2]. If this result is confirmed by another 
independent experiment, a vast unknown realm of the CDM cosmology will open up, 
because it is generally very contrived to build SUSY models producing excess positrons, 
but no excess antiprotons. It is clear from the above discussion though, that the decaying 
axino could be a possibility, since in that case the radiative or leptonic decay channels 
may be preferred over the hadronic ones El- In fact in the bilinear R-parity violating 
case, the two body decays into We~^ is either not open (for the hadronic axion models) or 
may be suppressed and therefore the antiproton flux from W fragmentations disappears 
leaving only the Z contribution. The direct channel into e~^e~i' may be dominant for 
models where the sneutrino v.e.v., i.e. the bilinear R-parity breaking, is suppressed and 
the trilinear R-parity violating couplings give the dominant decay. 

The model-independent case of a fermion decaying into e"'"e~i^ has been recently 
studied in the context of the PAMELA anomaly in [H], and there it was shown that 
such decaying DM particle may be a good fit to the data for lifetime around 10^^ s and 
mass above 300 GeV. The axino CDM could be a realization of this scenario with the 
appropriate choice of parameters. In that case though the axino has to be pretty heavy 
and therefore the reheating temperature very low. Even larger masses are probably 
needed for trying to accommodate the ATIC anomaly as well [15]. Note though that 
the PAMELA excess could be also due to astrophysical sources like pulsars |16] and 
then the PAMELA data give only a bound on the axino lifetime and R-parity breaking 
parameters. In that case also the possibility of a heavy axino > with a neutralino 
DM, which cannot fit very easily with the PAMELA anomaly, is not ruled out either. 



6. The LSPSM and colliders 



The signal of axino DM at colliders depends strongly on the nature of the LSPSM, 
which in turn depends on the SUSY breaking mechanism. In the constrained MSSM, 
where all the SUSY breaking parameters are derived by two common mass parameters, 
"^0, ^i/2i and a common trilinear coupling A at the unification scale, the value of tan j3 
and the sign of /i, the only allowed LSPSMs are the lightest neutralino and the stau. In 

+ The Icptophilic coupling for the axino was also observed in flipped SU(5) models [l5|[43]. 
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more general SUSY breaking models, of course other LSPSMs are allowed, in particular 
in the case of non universal Higgs scalar masses, the stop or the sneutrino. 

If the neutralino or the sneutrino are the LSPSM, it will be difficult to disentangle 
the two and prove that they are not DM. In both cases it would be necessary to measure 
their mass and couplings and realize that those parameters either give a too large DM 
energy density or are already excluded by direct DM searches. Then we would have 
good reasons to imply that the neutralino or sneutrino must be unstable on cosmological 
timescales, but it will be very difficult to determine what they are decaying into and if 
that includes the axino. Other, more indirect, collider signatures may arise in models 
with axino DM and 5*0(10) Yukawa unification [TF] . 

If the stau (or another charged sparticle) is the LSPSM instead, we will have the 
striking signal of an apparently stable charged heavy particle in the detector. In that 
case it will be clear that the LSP must be a very weakly interacting particle or that R- 
parity is violated, but we will need to measure and study the LSPSM decay to distinguish 
the two possibilities and identify if there is a DM candidate and which kind of particle 
it is. Unfortunately the astrophysical constraints on the R-parity violation scenarios 
discussed in the previous section seem to point to a quite long LSPSM lifetime, if the 
axinos are DM, and the decay would mostly happen outside the detector. 

6.1. How to distinguish the LSP from LSPSM decay ? 

The LSPSM decay can give information on the scenario and on the nature of the LSP, 
even if the LSP is not detected. In fact, the decay time and the branching ratios are 
model dependent and vary substantially e.g. between R-parity conserving and R-parity 
violating scenarios. In the first case, we expect that the dominant decay is the two- 
body channel into the LSPSM partner and the axino, while the next open channel 
the subleading radiative decay with an additional photon in the final state. If instead 
R-parity is violated, the LSPSM decays completely into SM particles with no missing 
energy apart for the light neutrinos. So for the case of a stau LSPSM, we have 

f ^ r a, r a 7 R-parity conserved; (22) 

T ^ r u^.fi Ur, b ^{h }fW~) R-parity violated; (23) 

therefore the R-parity violation case should clearly be visible via the large lepton number 
breaking, since e.g. the r and ^ final states arise from the same trilinear coupling and 
as well from the hadronic channel [H]. 

Moreover the angular distribution of the radiative decays into photon, a SM particle 
and missing energy, contains in general information on the spin of the LSP and the 
interaction vertex structure. This quantity can indeed play a key role in particular in 
order to distinguish between axino or gravitino LSP, that can give rise to similar NLSP 
lifetimes and similar "visible" decay channels [l9l [9]. In that particular case we will 
need to measure the branching ratio and the angular dependence of the radiative decay 
in order to reach a definitive identification [9]. 
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7. Conclusions 

We have discussed here different cosmological scenarios where the axinos play an 
important role in the DM question. If they are light, with masses in the MeV-GeV 
range, they can be the CDM if the reheating temperature is low and they are the LSP. 
In that case they can remain DM even if R-parity is broken, but the breaking has to 
be very suppressed. If instead axinos are heavy and not the LSP, they can still produce 
the necessary neutralino LSP abundance in their decay and dilute dangerous relics. 

In general, the presence of an axino LSP and DM relaxes many of the bounds on the 
super symmetric parameters, since the right number density of axinos can be obtained 
in a wider region of parameter space. Moreover the possibility of different LSPSMs and 
therefore very different collider signature arises. We expect LHC will soon clarify the 
situation. In the case of R-parity violation also astrophysical signatures could arise, but 
they are unfortunately strongly dependent on the axion model realization. 
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